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Abstract 

It has been shown that, for the two-qubit Heisenberg XY model, anisotropy and magnetic field 
may together be used to produce entanglement for any finite temperature by adjusting the external 
magnetic field beyond some finite critical strength. This interesting result arises from an analysis 
employing the Wootters concurrence, a computable measure of entanglement for two-qubit states. 
Recently, Mintert et al. proposed generalizations of Wootters concurrence for multipartite states. 
These MKB concurrences possess a mathematical property that enables one to understand the 
origin of this characteristic behavior. Here, we first study the effect of anisotropy and magnetic 
field on the multipartite thermal entanglement of a four-qubit Heisenberg XY chain using the MKB 
concurrences. We show that this model exhibits characteristic behavior similar to that of the two- 
qubit model. In addition, we show that this can again be understood using the same mathematical 
property. Next, we show that the six-qubit Heisenberg XY chain possesses properties necessary for 
it to have the characteristic behavior too. Most importantly, it is possible to directly measure the 
multipartite MKB concurrences of pure states. This may provide an experimental verification of 
our conjecture that for a Heisenberg XY chain of any even number of qubits, it is always possible 
to obtain non-zero genuine multipartite entanglement at any finite temperature by applying a 
sufficiently large magnetic field. 

PACS numbers: 03.67.Mn, 03.65.Ud 
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I. INTRODUCTION 



The one-dimensional Heisenberg models have been extensively studied in solid state 
physics (see references in jl|). Interest in these models has been revived lately by sev- 
eral proposals for realizing quantum computation [2] and information processing 3) using 
quantum dots (localized electron spins) as qubits. An intriguing phenomenon in such quan- 
tum systems with more than one component is entanglement. It refers to the non-classical 
correlations that exist among these components. Due to this nature, quantum entanglement 
has recently been recognized as an indispensable physical resource for performing classically 
impossible information processing tasks; such as quantum computation, teleportation |4|, 



key distribution [5| and superdense coding [6(. In any physically realistic consideration of 
quantum-information-processing device, it is important to take into account the fact that it 
would be in thermal contact with some heat bath. Consequently, entanglement in interacting 
Heisenberg spin systems at finite temperatures has been investigated by a number of authors 
(see, e.g., Ref.[?J and references therein). The state of a typical solid state system at thermal 
equilibrium (temperature T) is % = e _/3H /Z, where H is the Hamiltonian, Z = tre - ^ is the 
partition function, and (3 = 1/kT, where k is the Boltzmann's constant. The entanglement 
associated with the thermal state x is referred to as the thermal entanglement 

To achieve a quantitative understanding of the role of entanglement in the field of quan- 
tum information science, it is necessary to quantify the amount of entanglement that is 
associated with a given state. Several entanglement measures have been proposed. One 
famous example is the Wootters concurrence [8j. Consider a two-qubit state p, the Wootters 
concurrence 

C w [p] = max{Ai - A 2 - A 3 - A 4 , 0}, (1) 

where A& (k = 1,2,3,4) are the square roots of the eigenvalues in decreasing order of 
magnitude of the spin-flipped density-matrix operator R = p(a y '®a v ) p* (a y '®<j y ) , the asterisk 
indicates complex conjugation. Wootters concurrence is closely related to the entanglement 
of formation [9] . It appears as an auxiliary function that has to be evaluated when comput ing 



the entanglement of formation. Recently, Mintert, Kiis and Buchleitner (MKB) |10l . Ill ] 
proposed generalizations of Wootters concurrence for multipartite quantum systems, which 
can be evaluated efficiently for arbitrary mixed states. This is an extremely important 
development in our investigation of multipartite entanglement, a subject we have yet to 
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achieve a complete understanding. 

With the availability of such a good and computable measure of entanglement for systems 
of two qubits, many thoroughly analyzed the thermal entanglement in two-qubit Heisenberg 



models in terms of thermal concurrence 



12. 



13} 14, Q, Q, Q • The recent breakthroughs 



in the experimental physics of double quantum dot (see, e.g. , Ref.[l8j) shows that these 
studies are worthwhile pursuits. Interestingly, Kamta et al. [l5[] showed that the anisotropy 
and the magnetic field may together be used to control the extent of thermal concurrence 
in a two-qubit Heisenberg XY model and, especially, to produce entanglement for any 
finitely large T, by adjusting the external magnetic field beyond some finitely large critical 
strength. Such robustness is absent in the case of two-qubit Heisenberg XX model. A 
natural question is whether thermal entanglement of a Heisenberg XY chain with more 
qubits has similar behavior. This question is not only of practical importance but also of 
fundamental importance (see, for instance, Ref. 19| )• In this paper, we provide answers for 
Heisenberg XY chains with any even number of qubits by analyzing the MKB concurrences 
associated with their thermal states. 



Our paper is organized as follows. Since the results of MKB 10, LL1[ is crucial to our 
analysis, we will present them in the following section. However, we do so explicitly in the 
context of direct relevance to our paper, i.e., n-qubits. Specifically, when n = 2, the MKB 
concurrence coincides with the Wootters concurrence. It thus allows us to understand the 
special property of the thermal entanglement of the two-qubit Heisenberg XY model, in 
the light of a mathematical property of the MKB concurrences (see Eq.l fTTT) ). This together 
with a discussion of the two-qubit Heisenberg XY model [jjj] will be given in Section III. 
We show, in Section IV, that the thermal entanglement of the four-qubit Heisenberg XY 
chain has behavior similar to the two-qubit case. Namely, anisotropy and magnetic field 
may together be used to produce genuine multipartite entanglement [20|, |2l| for any finitely 
large T by adjusting the external magnetic field beyond some finitely large critical strength. 
This can again be understood in the light of the same mathematical property of MKB con- 
currences. Our results agree with those in Ref. 22|, which employed bipartite entanglement 
measure. In general, it is not sufficient to only consider bipartite entanglement. For instance, 
the Greenberger-Horne-Zeilinger (GHZ) state 23j is a state with genuine multipartite en- 
tanglement but yields zero entanglement between one particle and any other particle. On 
the other hand, the W state 24j is one where every particle is entangled with every other 
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particle, but it has no genuine multipartite entanglement [lCj, |2Cj. Multipartite entangle- 
ment surely is more interesting. Therefore, our study complements theirs. However, we 
must emphasize that our analysis more importantly demonstrates that the thermal entan- 
glement of the Heisenberg XY chains have the characteristic behavior due to two reasons. 
First, in the presence of a finitely large external magnetic field, the ground state of the 
model has non-zero genuine multipartite entanglement. Second, for non-zero temperatures, 
this ground state is the highest weight state (with weight ~ 1, i.e., the thermal state is 
almost a pure state) when a finitely large enough magnetic field is applied. In Section V, 
we show that this is indeed the case for the six-qubit Heisenberg XY chain. We thus have 
firm mathematical basis (Eq. fllip ) to establish that the thermal entanglement in this case 
will exhibit similar characteristic behavior. Most important of all, the MKB concurrences 
may be directly measured for multipartite pure states 
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26]. It is thus possible to experi- 



mentally verify our conjecture that for any even n, it is always possible to obtain non-zero 
genuine multipartite entanglement at any finite temperature by applying a sufficiently large 
magnetic field. Further discussions of this possibility and a summary of our results will be 
presented in the concluding Section VI. 



II. THE MKB CONCURRENCES 



Consider an n-qubit pure state \ip) G 'H\®'}ii ® ■ ■ •®7i n , the MKB concurrence is defined 
as the expectation value of a Hermitean operator A that acts on two copies of the state jlO ]: 

c[mm = vm®m^>®m- (2) 

In general, A could have the form 

A = E p^-^ij®^®--®^' (3) 
{-j 4 =±} + 

where p sl . S2 .... Sn . > 0, 

if = ^(n + + n+ + n -), 

p!P = ^n r (4) 

are projectors onto the symmetric and antisymmetric subspaces of Hj <8> TCj. Here, = 
(|00) ± |11))((00| ± (11|), nf = (|01) ± |10))((01| ± (10|) and {|0), |1)} is an orthonormal 
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basis of Tlj. The summation in Eq.Q is performed over the set {sj i = ±} + , which contains 
all n-long strings of +'s and — 's with even number of — 's. The superscript + indicates that 
the string with n +'s is not included in the sum. This is because for separable states, its 
expectation value in the symmetric twofold copy is non-zero. Terms with odd number of 
Pz s are naturally excluded in the sum since their expectation values in the twofold copy 
states is always zero. 

By choosing the value of all the p ai .a 2 .... Sn .'s in Eq.(j3]) to be 4, an entanglement monotone 
C n can be obtained. The resulting operator A n can equivalently be written as 4(7 — P^ <g> 



P^ <8> • • • <8> P+ ) [26j]. The concurrence C n of an n-qubit pure state \if>) can then be written 

as 

CnWW] = 2 1 - n/ y(2«-2)^|^)-^Trp J 2 . (5) 

The above summation runs over all (2 n — 2) reduced density operators pi of the state \ip). 
C n accounts for all possible types of entanglement in a state and takes the value zero if and 
only if the state is fully separable. 

For even number n of qubits, it is possible to define an MKB concurrence C (n) that detects 
multipartite entanglement, by choosing the operator A = = 2 n P^ <g> P^ <g> ■ ■ • <g> P^\ 
We note that when n = 2, 

c {2) u)(n = \m° v ®° y m = cwU)(n- (6) 

That is, the MKB concurrence coincides with the Wootters concurrence. And, for n — 4, 
we have 

C^[\tfj)(^j\] = | ^*\a y ® a y ® a y ® a y (7) 

which is also an entanglement monotone There is obviously no equivalent definition 
for the case of odd number of qubits since the expectation value of A^ are always zero for 
odd n. 

The MKB concurrence for a mixed state p of n qubits can be obtained via the convex 
roof construction: 

C\p] = inf j^PiC[|V>i)(V>i|], P = X>I^X^I j, (8) 

where the infimum is taken over all possible pure state decompositions of p. To evaluate 
C[p], consider the spectral decompositions p = ^ \4>i){4>i\ and A = \Xa)(Xa\, where the 



eigenstates \<fii) and \x a ) are subnormalized such that their norms squared are the eigenvalues 
corresponding to the states. If r is the rank of the operator A, it is possible to define r 
complex symmetric matrices T a with elements T? k = (<f)j \ ® (4>k\x a )- And, Eq.([8]) becomes 



C[p]=mf^Jj2\l VTaVT ^\ 2 ( 9 ) 

i V a 

where the infimum is now taken over the set of left unitary matrices V. 



It can be shown that the following inequality holds 



10|: 



C\p}>mfJ2\[VTV T U (10) 

i 

Here, the matrix r is defined to be Yl a z aT a i n terms of arbitrary complex numbers z a 
satisfying only the condition that Yl a \ z a\ 2 = 1- An algebraic solution of the inequality 
Eq. ffTUl) is given in Ref . |ly| to be max{0, Ai — J2j>i ^j} where Aj's are singular values of r 
written in decreasing order. For C^ n \ when A = A^ is of rank 1, T 1 = r and the lower 
bound in Eq. lflO]) turns out to be the exact value of C^ n \ 

In general, an optimization over z a is also necessary to obtain the optimal lower bound 
for C[p}. However, it is possible to obtain a good approximation to C[p] by approximating r 
with a matrix whose elements flO ] 

Ti i ~ rz 1 z ' \ ) 

V <^l| ® ® I0l) 

where |0i) is the eigenstate of p with the largest eigenvalue. This is the mathematical 
property of MKB concurrences that will play a critical role in our understanding of the 
characteristic behavior of the thermal entanglement of Heisenberg XY models. We will first 
illustrate this explicitly with the two-qubit Heisenberg XY chain in the next section. 



III. TWO-QUBIT HEISENBERG XY MODEL 

The Hamiltonian H2 for the anisotropic two-qubit Heisenberg XY model in an external 
magnetic field B m = rjj (77 is a real number) along the z axis is 

H 2 = + 7) Jo* ® of + ^(1 - j)Ja\ ® a\ + ^B m (af ® I 2 + h ® of), (12) 

where Jj is the identity matrix and crj, crj, cr| are the Pauli matrices at site j = 1,2. The 
parameter — 1 < 7 < 1 measures the anisotropy of the system and equals for the isotropic 
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XX model [l2] and ±1 for the Ising model [yj]. (1 + 7) J and (1 — 7) J are real coupling 
constants for the spin interaction. The model is said to be antiferromagnetic for J > 
and ferromagnetic for J < 0. The thermal concurrence associated with the thermal state 
X2, Eq. ffTSj) . can be derived from Eq. ffTTl) . It is invariant under the substitutions rj — > —rj, 
7 — > —7, and J — > —J. Therefore, we restrict our considerations to r? > 0, < 7 < 1, 
and J > 0. 

The eigenvalues and eigenkets of H2 are given by 



# 2 |$ u ) = £|$ u ), 

H 2 \^) = J]® 1 ), 

# 2 |$ 2 ) = -J|$ 2 ), 

# 2 |$ 3 ) = -£|$ 3 ), (13) 



where B = y-B^ + 7 2 J 2 = yrf + 7 2 J, 

l*°> = -==1^==[(B + S TO )|00>+7J|11>], 
l^ 1 ) = 7^[|01) + |10)], 
l$ 2 ) = ^[|01>-|10>], 

l$ 3 > = m _\- a W [(g " S m )|00> - 7 J|H>]. (14) 
The Wootters concurrence associated with the eigenkets, |$°) and |$ 3 ), are given by , 7 

\/V+7 2 

Hence, they represent entangled states when 7 7^ 0. We note that when rj = 0, |$°) and 
|$ 3 ) reduce to (|00) + \ll))/y/2 and (|00) - |ll))/\/2 respectively, so that the eigenstates 
are the four maximally entangled Bell states: I^Beii)) I^Beii); I^Beii)) an( l I^Bcii)- And, in 
the limit of large 77, 

C w [\^)(^\]=C w [\^)(^\]^lV-\ (15) 

only going to zero asympotically when rj is infinitely large. In contrast, when 7 = 0, 
|$°) = 1 00) and |$ 3 ) = |11) are product states with eigenvalues r/J and —77 J respectively, 



though I*!? 1 ) and |$ 2 ) remain the same 12] 



For the above system in thermal equilibrium at temperature T, its state is described by 
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the density operator 

3 

X2 = 52wi\&)(&\ 

i=0 

= L[ e -W\$0)($°l + e - /3J |$ 1 )($ 1 | + e /3J |$ 2 >($ 2 | + e /3S |$ 3 >($ 3 |], (16) 
z 2 

where the partition function Z2 = 2cosh/3£> + 2 cosh j3 J, the Boltzmann's constant k = 1 
from hereon, and /3 = 1/T. After some straightforward algebra, we obtain 

Ax = 

^2 

A2 = ic"^, 
^2 



1 

^2 



1 + 



2 7 2 J 2 
£ 2 


sinh 2 /3£ + 


2 7 J 
B \ 


A + 781/8 

/ B 2 


sinh 2 /3Bsinh/3B, 




2 7 2 J 2 
£ 2 


sinh 2 (3B - 


2 7 J 
B \ 


A + 7272 

/ B 2 


sinh 2 /3Bsinh/3B. 


(17) 



A 4 = ^1/1 + 

In the zero-temperature limit, i.e., /3 — > 00, at which the system is in its ground state, 
Eq. (FUJI) reduces to the following three possibilities. 



(a) < r] < y/1 -7 2 : 

X2 = _l[ e /"|$ 2 )($ 2 | + e ^|$3)($3|] _> |$ 2 )($ 2 |, (18) 
Z2 

with 2" 2 = e /3J + e /3S . Equation (fTTl) gives CwlXz] = 1> its maximum value, in agreement 
with the fact that |$ 2 ) is a maximally entangled Bell state. 

(b) rj = v/T^y : 

X2 ^i[|$ 2 )($ 2 | + |$ 3 )($ 3 |]. (19) 
From Eq. flTTI) . the above equally weighted mixture has 

C w [ X 2] = \{l-l). (20) 

(c) r] > a/1 - 7 2 : 

X2 ^|$3 )($ 3| ; (21) 

and Eq. fTlT|) yields accordingly 

C w ix 2 ] = J— . (22) 
V V + T 
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Therefore, for a given 7, r7 critica i = a/1 — 7 2 marks the point of quantum phase transition 
(phase transition taking place at zero temperature due to variation of interaction terms in 
the HamUtcuan of a 8y8 te,„ fl). For va.oe S of 7 other than 7 = i, there is a 8 „dde,t increase 
or decrease in C[x2\ at 77 cr i t i ca i, depending on whether 7 > | or 7 < |, before decreasing to 
zero asymptotically, as r\ is increased beyond the critical value r/ critical [15|. Here, we focus 
on the behavior of the model at non-zero temperatures and subject to magnetic field of 
appropriate strengths. 

At non-zero temperatures, due to mixing, C\y[X2] decreases to zero as the temperature T 
is increased beyond some critical value. In fact, as T — > 00, the statistical weights Wi — > 1/4 
for all i and Cw[X2] - > 0. However, for a large but finite T, 

e& = 1 

can always be made as close to unity as possible by increasing the strength 77 of the external 
magnetic field. That is, when 77 is large enough, only |<3> 3 ) contributes significantly to the 
thermal state X2- The 77 required for this to occur depends on T, larger rj for higher T. This 
is always possible for finite T. 

Next, we note that in the limit of large 77, 

A 1 ^A 2 ^0^A 4 , X 3 ^jr ] - 1 . (24) 

It follows that 

C w [ X2 ]«7^ 1 ~C w [|$ 3 )($ 3 i]. (25) 

Hence, for a finitely large T, the thermal concurrence of the system in a large enough 
magnetic field is very well approximated by the concurrence of |$ 3 ) in the same magnetic 
field. In other words, the entanglement associated with the thermal state X2 m this case is 
mainly due to that associated with the eigenstate |$ 3 ). 

In summary, for any non-zero T and an appropriate rj, |$ 3 ) is the highest weight eigen- 
state. It follows from Eq. fllll) that this is the state which will mainly contribute to the MKB 
concurrence C^ 2 - ) [x2] or equivalently the Wootters concurrence (see Eq.([n])). An important 
point to note here is that |$ 3 ) has non-zero concurrence as long as 77 is finite (see Eq.f ll5l) ). 
Another is that 7 7^ 0, otherwise |$ 3 ) will be a product state with no entanglement. These 
clearly explain the characteristic robustness of the thermal entanglement of the two-qubit 



Heisenberg XY model described in Ref 



15j - the system at finite T can always be entangled 



provided large enough magnetic field is applied. They also identify the necessary character- 
istic features for a model to exhibit such behavior. We shall illustrate that this is indeed 
the case for the four-qubit Heisenberg XY chain in the next section. 



IV. FOUR-QUBIT HEISENBERG XY MODEL 

The Hamiltonian H n for an anisotropic n-qubit Heisenberg XY chain in an external 
magnetic field B rn = r)J along the z-axis is 

T n 

#n = 3 £[(* + ^H+i + (1 " 7)^ + i + m-i (26) 
i=i 

where the periodic boundary condition cr^+i — °\ ( a — x -> Vi z ) applies. Like in the two-qubit 
model, we consider the case when 77 > 0, < 7 < 1, and J > 0. 

In this section, we consider n — 4. After some straightforward algebra, we obtain the 
eigenvalues and eigenvectors of H±. Firstly, we present 



where 



and 



CO 



# 4 |$ 15 ) = -cj+J|$ 15 ), (27) 



± = y/2 \J [ V 2 + 2(1 + 7 2 )] ± v^f + 2(1 + 7 2 )] 2 - 8^, (28) 



|$ 15 ) = AT-(fi-|0000) + ^ |0011) + ^3 |0101) + ^2 I0110) 

+Q2 |1001> + Q3 |1010) + or |noo) + |1111>), (29) 

With n± = i/y^TItW+WW+WW?, 

(2?7 ± uj+){uj +2 - 8) - 87 2 (?? ± u + ) 



87 2 r7 
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Slf = ± 2 V±f . (30) 

|<3>°), which can be obtained from |<3> 15 ) by substituting f2~ {% = 1, 2, 3) with fit or fl~ — > f2+, 
satisfies i?4|<E>°) = c^Jl^ ). The MKB concurrences, as defined in Section II, for |$ 15 ) are 
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given by 

C (4) [|$ 15 )($ 15 |] = 2 fi r + 2(^ 2 ) 2 + (^) 2 



C 4 [|$ 15 )($ 15 ' 



i + (fiD 2 + 4(fi 2 -) 2 + 2(fij) 2 ' 

V7[2(fi 2 -) 2 + (J^) 2 ] 2 + 2[4(^) 2 - + 4] [2(fi 2 ) 2 + (ftj) 2 ] + 7(«D 2 



i + (fir) 2 + 4(fi 2 ") 2 + 2(n 3 -) 2 



(31) 



In the limit of large f], 



+ ~ 9 2 7 2 4 7 2 - 7 4 
u; ~ 2?7 H h 



77 r/ 3 
7 2 2 7 2 - 7 4 



2r/ 2 2?7 



4 ' 



n 2 « + 

2 2?7 4?7 3 ' 

7_ 4 7 - 3 7 3 
r/ 2 2r/ 4 



It follows that 



C (4) [l$ 15 )($ 15 |]-^ + ^i^, (33) 



C 4 [|$ 15 )($ 15 |]^^ + 24 ^ ^\ (34) 
77 Arf . — . 

We note that |$ 15 ) is a state with genuine four-partite entanglement [27j, which remains 



non-zero even for large 77 and going to zero only in the asymptotic limit of infinite magnetic 
field. 



Secondly, we have 



where 



and 



# 4 |$ 14 ) = _[(«+ + a -) 7 + 2 ]J|$ 14 >, 

tf 4 |$ 13 > = -[(« + + cT)7 - 2] J|$ 13 ), (35) 



a* Vg+g±^ ; (36) 



|$ 14 ) = — (-q~|0001) + orlQOlO) - q~|0100) + 10111) 

2^1 + (a-) 2 

+oT|1000) - |1011) + |1101) - |1110>), 
|$ 13 ) = ^ = = (-a"10001)-a"|0010)-g~|0100) + 10111) 

-aT|1000) + |1011) + 1 1101) + 1 1110>). (37) 

11 



Corresponding to these states are and |$ 2 ), which can be derived from |$ 14 ) and |$ 13 ) 
respectively by -or -> a + . They satisfy # 4 |$ 1,2 ) = [(« + + a~h ± 2] J|$ 1,2 ). The MKB 
concurrences for |$ 14 ) are given by 

C ( 4 )[|$ 14 )($ 14 |] 2n ~ 



1 + (ar 



CHO^I] = _J_; 3 + 8( a -)^3(a.-)^ (3g) 



Hence, |$ 14 ) is also a state with genuine four-partite entanglement [20J, |2l|. It reduces to a 
W state when 7 = 0. 

Thirdly, we have # 4 |$ 3,4 ) = r?J|$ 3 ' 4 ) and H A \^ n > 12 ) = -rjJ\<$> n > 12 ), where 

|$ 3 ' 4 ) = — (|0001) ± |0010) — |0100) =F 1 1000)), 

|$ n ' 12 ) = — (|0111) ± 1 101 1) — 1 1 101) =F 1 1 1 10) ) . (39) 

These states belong to the family of W states which are entangled but do not contain genuine 
four-partite entanglement 20] . Fourthly, by substituting uj + in Eq. flBUl) with u~ we obtain 
the corresponding , A^ and A^ in terms of which we express 

|$io,ii) = jV±(A± |0000) + Af |0011) + Af|0101) + Af |0110) 

+A±|1001) + Af |1010) + A^|1100) + 1 1111)). (40) 



Here, N% = 1/^1 + (Af ) 2 + 4(A±) 2 + 2(A±) 2 . They satisfy tf^ 10 - 11 ) = ±aT J|$ 10 ' n ). 
Lastly, we have the following four degenerate eigenstates with eigenvalue zero: 



$ 6 ) = 


^(10011) 


-11100)), 




$ 7 ) = 


^(10101) 


-11010)), 




$ 8 ) = 


^(loiio) 


-11001)), 




$ 9 ) = 


1 (10011) - 


|0110) - |1001) + |1100>) . 


(41) 



The eigenstates |$ 6 ' 7,8 ) belong to the family of GHZ states that contain only genuine four- 
partite entanglement, while |$ 9 ) is a product state of two Bell states. 

As in the case of the two-qubit model, we construct the thermal state of the four-qubit 
Heisenberg XY chain as follows: 

15 

XA = Y,w l \V){V\ = -e-P H \ (42) 
12 



Here, the partition fuction 



Z 4 = 4 + 4 cosh (3t]J + 2 cosh (3[(a + + a~)^ + 2] J 

+2 cosh p[(a+ + a~)7 - 2] J + 2 cosh ^ + J + 2 cosh /3a; ~ J. (43) 

At non-zero temperatures the state of the system becomes a mixture of the energy eigenstates 
with statistical weights 



3 /3[(a++a-) 7 +2]J 



«?14 



(44) 



A. Zero Temperature 



In the /? — > oo limit, 



X 4 = — [ e ' 3 [( a++Q "^+ 2 ] J |$ 14 )($ 14 | + e ^ +J |$ 15 )($ 15 |], 



(45) 



with Z A = e^ a++a ) 7+2 l J + J . Like in the 2-qubit Heisenberg XY model, there are thus 
two possible lowest energy states, namely |$ 14 ) or |$ 15 ), depending on the strength of the 
applied magnetic field. As rj is increased from zero, there are in general two instances when 
(a + + a~)7 + 2 = u + . We let 771 and r\ 2 denote the solutions. Their dependence on the 
anisotropy of the system are plotted in Fig. HJ As 7 is increased from zero, both r\\ and 772 
become smaller and converge to zero when 7 = 1. In fact, for the Ising model (7 = 1), both 
|$ 14 ) and |$ 15 ) are the lowest energy states of the system. The system state is an equal 
mixture of both states: 

X4 = ^[|<f 14 )($ 14 | + |$ 15 )($ 15 |]. (46) 

But, once the external magnetic field is turned on, |<3> 15 ) becomes the only ground state of 
the model regardless of the strength of the field. 

It follows that for < 7 < 1, depending on 77 we have the density operator of the system 



X4 



|$ 15 )($ 15 | 
1[|$15 )( $15| 

|$ 14 )($ 14 | 
1[|$14)($14| 

|$ 15 )($ 15 | 



< 7] < Tji 
|$ 14 )($ 14 |] T) = 771 

Vi < V < V2 

|$15)($15|] v = r]2 
7] > 7] 2 



(47) 
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FIG. 1: The two transition 77's for the system at zero temperature plotted against the anisotropy 
7- 

We can therefore calculate the MKB concurrrences of the system at zero temperature {(3 — > 
00) for the different regions of magnetic field strength. Plots of and C4 against magnetic 
field for different values of anisotropy 7 are shown in Fig. [2] and [31 For both concurrences, 
there are sharp changes at the transition magnetic fields due to quantum phase transition. 



C !4) 




• ' — 1 1 1 1 — 1 1 

0.5 1 1.5 2 2.5 3 



FIG. 2: The MKB concurrence which is an entanglement monotone, plotted against magnetic 
field i] for different values of anisotropy 7. In general, there are sharp changes in at the two 
transition fields rji and 772- 

B. Non-zero Temperatures 

In general, as the temperature of a system is increased, its density operator becomes 
closer to the maximally mixed state, -I, where n is the dimension of the Hilbert space 
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FIG. 3: The MKB concurrence C4 plotted against 77 for different values of 7. Similar to Fig|2j there 
are sharp changes of the concurrence at the two transition magnetic fields. 

and / is the identity operator. There thus exist critical temperatures T c beyond which 
the MKB concurrences of the system become zero, like in the two-qubit case (see Fig. H]). 
The existence of T c 's is guaranteed by the fact that a state becomes separable when it is 
sufficiently close to \l [281 ]. 
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FIG. 4: The points where the MKB concurrence equals zero are plotted for two different values 
of 7. The region above each curve is the region where = 0. 

On the other hand, it can also be seen from Fig. H] that it is always possible to have 
non-zero entanglement in the system by applying a sufficiently large magnetic field to it. 

15 




x 



In particular, even when T > T c , one could reintroduce entanglement into the system by 
increasing rj. This is always possible for finite temperatures as long as 7 ^ 0. Hence, the 
thermal entanglement associated with the four-qubit Heisenberg XY model exhibits the 
same characteristic robustness as in the two-qubit case. We note that, in contrast to the 
two-qubit case (see Fig. 3 in Ref.[15|), each of our graphs in Fig. H] has two "singular 
turning points". This is due to the fact that there are two transition 77's instead of one in 
the two-qubit model. 

In order to understand this robustness, we draw inspiration from the two-qubit case. We 
observe that for any finite temperature T, the statistical weight W15 when 77 is large enough, 
is given by 

w 15 = — exp (3uj + J = \ (48) 
z 4 ^ 

where £ = 1 + e~ 2(3uJ+ 3 + Aq~^ +j + 4 e -^ +J cosher/ J + 2e~^ +J cosh/3[(a+ + a~)j + 2} J + 
2e~ l3uJ+J cosh/5[(a + + af~)7 — 2]J + 2e~ l3u)+ 3 cosh/3u; _ J. It can be shown that the large 77 
behavior of W15 is independent of 7. Furthermore, by increasing 77 appropriately, £ can be 
made to be as close as possible to unity. This results in |$ 15 ) being the only state that 
significantly contributes to the thermal state x<i- m the light of Eq. fflT]) . we may conclude 
that 

C (4) [x 4 ] ~C (4) [|$ 15 )($ 15 |]. (49) 

That this is indeed the case has been established numerically (see Table [I]) . Since the MKB 
concurrence can be determined exactly, it is calculated for states %4 an d |$ 15 ) under 
different combinations of 77 and T (some results are shown in Table [I]) . The two values for 
each combination of rj and T are then compared. For a given temperature, the two values, 
C^[x4] and C^ 4 - ) [| < I )15 )( < I )15 |], agree when large enough magnetic field is applied. 

Hence, for 7 7^ 0, the revival and robustness of entanglement in the Heiseberg four-qubit 
model can be understood, as in the two-qubit case, in terms of the large 77 behaviors of both 
the MKB concurrence C (4) [|$ 15 )($ 15 |] (Eq.(|I5J) and the statistical weight w 15 (Eq.pBl). We 
may conclude, when u> 15 w 1, that the entanglement associated with xa is of the genuine four- 



partite kind 
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2l|. The total entanglement as measured by C 4 thus also undergoes revival 



since genuine four-partite entanglement is only one kind of the entanglement measured by 
C4. In contrast, for 7 = and large rj, the ground state is the product state 1 1111) . This is an 
important distinction between the four-qubit XX model and the XY model. Consequently, 
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|$ 15 )($ 15 | 


1 





1 


0.0000180069 


0.0000180069 


1.79177xl0~ 7 


1.79177xl0" 7 


5 





1 


0.0000180068 


0.0000180069 


1.79177xl0~ 7 


1.79177xl0" 7 


10 





1 


0.0000174316 


0.0000180069 


1.79177xl0~ 7 


1.79177xl0" 7 


50 





1 





0.0000180069 


1.79175xl0~ 7 


1.79177xl0" 7 


100 





1 





0.0000180069 


1.07513xl0~ 7 


1.79177xl0" 7 



TABLE I: A comparison between C^[ X a] and C( 4 )[|$ 15 >($ 15 |] for some combinations of rj and T. 

no revival of entanglement is observed for the case of four-qubit isotropic Heisenberg XX 
chain. 

We may apply the above analysis employing C^ n ' to study Heisenberg XY chains with 
even number n of qubits. Firstly, we determine if the ground state remains genuinely 
multipartite entangled at large r]. Secondly, we determine if the ground state statistical 
weight w g can be made very close to 1 at large rj. Through Eq. liTTj) . these two properties 
together imply that for any finite temperature, nonzero genuine multipartite entanglement 
can always be obtained by applying a large enough magnetic field. We show that the six- 
qubit Heisenberg XY chain has both properties in the next section. 

V. SIX-QUBIT HEISENBERG XY MODEL AND BEYOND 

In this section, we show that the six-qubit Heisenberg XY chain (with 7^0) does 
indeed possess the characteristic properties, which enable the model to have non-zero thermal 
entanglement at any given finite temperature when subject to an external magnetic field of 
appropriate strength. To this end, we determine the eigenstate of H 6 whose eigenvalue 
Eg = —XJ is the minimum when 7] is large. Here, 

A = ^3(2 + 2 7 2 + ri 2 ) +2k + 2^ V W + ^ 2 )((3 + 7 2 + V 2 ) + «0 (50) 
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with k = V7 4 + (-3 + f) 2 ) 2 + 2 7 2 (3 + rj 2 ). And, 



= ^(1000000)0! + |OOOO11)0 2 + |OOO1O1)0 3 + |OOO11O)0 4 + |OO1OO1)0 4 + |OO1O1O)0 3 



with 



and 



+ 
+ 
+ 
+ 
+ 



|OO11OO)0 2 + |OO1111)0 5 + |O1OOO1)0 3 + |O1OO1O)0 4 + |O1O1OO)0 3 + |O1O111)0 6 
|011000)©2 + |O11O11)0 7 + |O111O1)0 6 + |O1111O)0 5 + |1OOOO1)0 2 + |1OOO1O)0 3 
|1OO1OO)0 4 + |1OO111)0 5 + |1O1OOO)0 3 + |1O1O11)0 6 + |1O11O1)0 7 + |1O111O)0 6 
|110000)©2 + |11OO11)0 5 + |11O1O1)0 6 + |11O11O)0 7 + |111OO1)0 5 + |111O1O)0 6 
|1111OO)0 5 + |1111H)08) (51) 



0i 

©2 

©3 

©4 

©5 

©6 

7 
©8 

N 



{24J 2 - (A - 3Jr]){\ - Jri)}®8 + 2J 2 (A - Jr])( + 8J 2 r 



2 J 7 (A + 3 J77) 



(A + 37^)0! 
6J7 

(A + 3Jr?)0 8 + (A - Jr])r + 2J 2 ( 
6 J7 2 

3(2 - 7 2 )0 8 + (A - Jr])( + 2r 
3 7 2 
(A -3Jt7)6 8 



6J7 
30 8 + r 



~ 3 7 ' 
(A - 3J??)08 - 2J 2 C 
6J7 ' 

A 4 - 2 J 2 A 2 (2 + 2 7 2 + rf) + J 4 {V (?? 2 - 12) + 4 7 V + 4)}, 

1 

V© 2 + 2 + 3(20 2 + 20 2 + 2 + 20 2 + 20 2 + 2 ) ' 



T = -8J\r]{\ 2 - J 2 (6 - 2 7 2 + r7 2 )}, 

( = 4A{A 2 ( 7 2 - 1) - J 2 (4 7 4 - 9 V 2 - 4 7 2 + 7 V)}. 

The MKB concurrence C^ 1 is calculated for this state, giving 

C (6) [|$ s )($ s |] = 2iV 2 |0 1 8 + 6020s + 60306 + 30 4 7 | 



In the limit of large 77, 



(52) 



-3 



(53) 
(54) 



(55) 



(56) 
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It can again be shown that is a genuine six-partite entangled state [20l. l2lj. In addition, 
the statistical weight w g of the ground state \Q 9 ) can be shown numerically to be close 
to unity when an appropriately large magnetic field is applied. The ground state \§ 9 ) 
therefore possesses the two desired properties for the system to have non-zero genuine six- 
partite entanglement at any finite temperature, provided an appropriate magnetic field is 
applied. We conjecture that the robustness of genuine multipartite entanglement is a general 
property of Heisenberg XY chain with even number of particles. This conjecture can be 
experimentally tested, as will be discussed in the next section. 



VI. CONCLUSIONS 



In this study, we have investigated in detail the origin of the robustness of genuine mul- 
tipartite entanglement in two-, four-, and six-qubit Heisenberg XY models. Two important 
properties possessed by the ground states |$ 9 ) of these models, which enable them to ex- 
hibit the characteristic robustness, were identified. The first property, namely the statistical 
weight w g associated with can be made very close to unity by applying a large enough 
magnetic field, allows us to use Eq. fllll) and conclude that the MKB concurrence of the 
thermal state \ equals that of the ground state. It follows that we could obtain fairly ac- 
curately the MKB concurrence of \ (a mixed state) by calculating the MKB concurrence 
of (a pure state) in this case. The second property being that remains genuinely 
multipartite entangled under such a magentic field, then guarantees that there is non-zero 
genuine multipartite entanglement for any finite temperature as long as sufficiently large 
magnetic field is applied. These properties allow us to extend our study to XY chains with 
any even number of qubits. Heisenberg XY chains with odd number of qubits were not 
considered in this study because = for odd n. However, one could similarly study the 
robustness of other kinds of entanglement employing the other MKB concurrences. 

In fact, the applicability of our analysis is not only restricted to the Heisenberg XY 
chains. As long as the ground state of a system can be made to dominate the thermal state 
at any finite temperature by adjusting some parameters of the system Hamiltonian, Eq. ffTTj) 
tells us that the MKB concurrence of the thermal state is given by that of the ground state. 
So, if in addition, the ground state remains multipartite entangled under these conditions, 
similar robustness of multipartite entanglement is expected to be observed. Therefore, our 
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analysis can be used to identify possible candidates for realization of quantum computation 
at finite temperatures. 

The MKB concurrences C^ n > and (L can be directly measured for pure states, if two 



copies of the states are available 
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261 ] . This is obviously due to the fact that the MKB 



concurrences are defined in terms of expectation values of Hermitean operators. Indeed, 
this has lead to the direct measurement of Wootters concurrence, previously thought to be 



not directly measurable, for pure states in laboratory [25[. Since we are interested in the 
region where the ground state |<3> 9 ) is the only state that contributes significantly to the 
thermal state Xi the state x is "almost pure" and therefore it is possible to measure the 
MKB concurrences of x with a high degree of success. The results obtained here can thus 
be experimentally verified. 

In conclusion, we have established a rather general method of identifying systems that 
exhibit robustness of multipartite entanglement at finite temperatures. Our analysis rest on 
the mathematical property of the MKB concurrences, namely Eq.t flTT) . The MKB concur- 
rences can be directly measured and therefore experiments can be carried out to verify the 
results of any analysis. 
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